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Abstract 

Cooldown filni-boiling-heat transfer data from small diameter spheres (0.0397 to 
1.27 cm) are presented, and the effects of curvature and similarities of the theo- 
retical analyses for the flat plate, cylinder, and sphere are assessed. The data 
are in good agreement with theory but ai’e significantly higher than predicted from 
the classic Bromley equation. 

INTRODUCTION 

Film-boiling heat transfer from small diameter spheres-is of fundamental interest 
in studying the effects of curvature, and of practical importance in engineering and 
biology. In nuclear reactors, small fission fragments (approximated by a sphere) 
can get "loose'* and film boil in the coolant passages producing voids. In preser- 
vation of biological specimens and human organisms, the heat transfer rates of 
snip 11, fluid globules must be accurately known to control the freezing process. 

Film boiling heat transfer to a flat plate or a surface with no curvature has been 
studied by several authors (e.g. , 1,2), and the fundamental relation presented by 
Berenson (Ref. 1) 


Nu £ = 0.424(Rap 1 / 4 

correlated much of the available data. Film boiling from a cylinder (single curved 
surface) has also been studied by several authors (e.g. , 3,4). Extending the pre- 
vious work, Breen and Westwater (Ref. 4) empirically correlated much of the avail- 
able cylinder data 


Nuj = (0.069\ c /d + 0.59)(Ra|f A c ) 1/4 

Compared with the classic Bromley correlation, (nUcI - 0. 62(Racj)* // ^], ref. 3, this 
correlation predicts a significant increase in heat transfer for very small and very 
large diameter cylinders. Later, Baumeister and Hamill. (Ref. 5) analyzed film 
boiling from the flat plate and the cylinder and were able to correlate available data. 
Their model extended the understanding of the film boiling process. 

Film boiling from a sphere has been studied by references 6, 7, 8, 9, and 10. The 
analysis of Frederking and Clark (Ref. 10) correlates much of the available data for 
large diameter spheres. 


Nu d « 0. 14 Ra£ l /3 


Recently, Hendricks and Baumeistcr (Ref. 7) analyzed film boiling from both small 
and large spheres and were also able to correlate these data. However, a signifi- 
cant difference exists in the correlations for the small-sphere case. Thus the pur- 


pose of this paper is to present experimental film-boiling heat transfer data for small 
spheres, test the analysis of reference 7, and also to assess the effects of curvature 
and similarities of the theoretical analysis for the flat plate, cylinder, and sphere. 
The analyses are outlined in the Appendix and the results given in table I. 

APP ARATUS AND PROCEDURE 

Cool down heat transfer data were taken on six different spheres with diameters 
ranging from 0.0397 to 1.27 cm. In each experiment, room temperature spheres 
were quenched in a bath of liquid nitrogen. 

The smallest spheres (0.0397 and 0. 0794 cm) were Tungsten-Carbide while the four 
remaining larger spheres (0.1585; 0.318; 0.397; and 1.27 cm) were copper. Tung- 
sten Carbide was the only commerically available material for the small diameter 
spheres. The 0.0397, 0.0794, 0.1585, and 0.397 cm spheres were instrumented 
with 0.00126 cm diameter Chromel-Alumel thermocouple; while the 0.318 and 
1.27 cm spheres had none. The heat transfer to the uninstrumented spheres was 
evaluated by measuring the quenching time (discussed in detail later). The couples 
were spot welded to the Tungsten -Carbide spheres and embedded into the copper 
spheres. The thermocouple output was recorded and reduced; a typical plot is shown 
in Figure 1. The accuracy of the recording system was within ±3° K and the repro- 
ducibility of the data was within ±10%. 

In Figure 1, tL e ^, represents the time for film boiling to cease, while tQ, repre- 
sents the time required for the sphere to reach the bath temperature. A sharp break 

in the time temperature curve can be seen in Figure 1 at and condition (1) 

follows: 


'■Leid w (1) 

The measurement of Iq will be euown to give a good first order estimate of the heat 
transfer coefficient. For the noninstrumented spheres, the quenching time was mea- 
sured using a 1/100 second stopclock and coordinated visually. Each sphere was 
quenched ten times and an average time was taken. The deviation ranged from ±15% 
for the 0.0397 cm sphere to less than 1% for the 1. 27 cm diameter sphere. 

THERMOMETRY ERROR 


The thermocouple conduction errors in measured surface temperature were esti- 
mated by modifying the analysis of reference 11, and found to be less than 1° K; 
no correction was applied to the data. 

The sensible heat loss by conduction through the thermocouple wires was estimated 
using the cylindrical fin equation 

“TC “ ! "AV u < T w - V - 2 ™ u ,^ e k L N “L < T w - V <2> 

In these calculations, values of k e == 0. 044 cal/cm-sec-K and Nu L - 1 were as- 
sumed and q,p£ altered the heat transfer coefficient as much as 25% for the 0.0397 
cm sphere. 




EXPERIMENTAL HKAT TRANSFER COEFFICIENT 


The energy equation for quenching' a sphere is given as 


-peVD t (T w - T s ) hA(T w - T s ) + q TC (3) 

whom the bath is at the saturation liquid temperature T s . At any instant the thermal 
capacitance (pcV) and the wall temperature T w are assumed uniform throughout the 
sphere. The latter assumption is justified because the value of the F 4 modulus 
(Bj - hd/k) is low. By rearranging equation (3), the film boiling heat transfer coeffi- 
cient can be ex pressed as 

h -{pcV/A)D t ln(T w - T g ) - q Tr /A(T w - T fl ) (4) 

Since the geometry and properties are known, the heat transfer coefficients we”e de- 
termined directly from equation (4) using values of ln(T w - T fl ) from the experi- 
mental data of Figure 1. The results are illustrated in Figure 2, For the 0.0397 and 
Or 07 97 cm spheres two heat transfer paths were found. The reason for these paths 
is-at present unknown, and attributed to uncertainities in the data. 

Average Heat Transfer Coefficient 

As seen in Figure 2, for AT > 90° K the heat transfer coefficients are reasonably 
constant. Consequently, assuming average properties and neglecting heat losses, 
an average heat transfer-coefficient is found by solving equation (3). This gives 

< T W " T S> /T W 0 ■ T s> ' exPH/V “ ex P(~thA/pcV) (5)- 

where T Wq is the initial sphere temperature and r c is the sphere time constant. 

For a given thermal resistance (l/hA)and capacity (pcV), the system time constant 
r c can be specified. To determine r c , the average heat transfer coefficient, F, 
was estimated from the data in reference S. Assuming condition (1), the quench- 
time to the calculated-time-constant ratio tQ/V c was experimentally found to be 
1.85; consequently a first-order estimate, of the average heat-transfer coefficient 
can be made 


h y - (pcd/e)/(t Q /l . 85) (6) 

Without any attached instrumentation, simply observing the quenching time, equa- 
tion (6) gave values for the heat transfer coefficient which were in good agreement 
with the values determined from equation (4). 

DISCUSSION OF THEORETICAL AND EXPERIMENTAL RESULTS 

In the derivation of the heat-transfer coefficients outlined in the Appendix, the heat 
transfer coefficients (listed in table Ij-arc found to be a function of two key similarity 
parameters f and Bond number B Q . The theoretical expressions for Nusselt num- 
ber are plotted as a function of Bond number in Figure 3 for the flat plate, cylinder, 
and sphere along with the classical Bromley equation. As is evident, at low Bond 
number a significant enhancement in heat transfer due to curvature is predicted for 


both the cylinder-and the sphere, lb ing these same parameters, the experimental 
heat-transfer data for spheres aro compared with theory. In. Figure 4. 

It is apparent that the data significantly depart from the large sphere analysis of _ 
reference 10 as woll as the classic Bromley relation and follow those trends estab- 
lished in Figure 3. However, while the analysis of reference 7 predicts the gener- 
al data trend, the 0.0397 cm sphere data appear on the average to be about 25% 
higher than theory. Three factors which may contribute to this deviation are: 

(1) difficulties in assessing large time-temperature gradients, (2) significant heat 
loss through t' c thermocouple, and (3) uncertainty in the heat transfer data as il- 
lustrated in Figure 2 and the specific heat of tungsten carbide. The vertical scat- 
ter in the data of Figure 3 exhibits a dependence on temperature difference which 
is not accounted for by the theory. (This is also true for the flat plate and the 
cylinder.; Perhaps a variable property analysis would encompass these trends. 

In addition, data for a 10.16 cm diameter sphere (Ref. 8) and those collected in 
refei’ence 7 are in agreement with the theory as seen in Figure 4. These data verify 
the~theory for a Bond number range from 0. 132 to 8700. 

SIMILARITY 


Similarities exist between the flat plate, cylinder, and sphere in the method of 
modeling, governing equations, boundary conditions, constraints, and evolved 
similarity parameters. Sufficiently allied are these criteria as to suggest they 
form a solution basis for all pool film boiling problems (see Appendix for details). 

Some-physical similarities between the sphere and cylinder are noted in Figure 5. 
At small B 0 , the vapor release is periodic and the field quiescent. The transition 
regimes are similar and note the analogous chaotic appearance at large B Q . Fig- 
ure 6 shows the similar idealized vapor dome configurations chosen to model the 
phenomena. The vapor dome spacing pa rameter \ is related through hydrody- 
namic stability theory to B 0 . 

CONCLUSIONS 


(1) The theoretical analysis of reference 7 predicted a significant increase in the 
film -boiling heat -transfer coefficient for small-diameter spheres, compared with 
previous analysis for large-diameter spheres. Experimental film boiling data for 
small spheres ranging from 0.0397 to 1.27 -cm confirm the theoretical prediction. 

(2) The Bond number is shown to be a key correlating parameter in relating the en- 
hancing effect of system curvature on the film boiling heat transfer coefficient. This 
enhancement becomes extremely pronounced at small Bond numbers (small spheres 
or cylinders) and approaches that of the flat plate at large Bond numbers. Fquiva- 
lently, if the radii of curvature are small, or very large, compared with the vapor 
dome spacing parameter £ , the effects of curvature augment the heat transfer over 
that predicted by the Bromley relation. 
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APPENDIX A 

General Method of Solution 

The method of modeling and the theoretical solutions for film boiling from the flat 
plate, cylinder, and sphere are quite similar. In this Appendix we shall review the 
models, the governing equations, their solutions and seek similarity forms which 
will enable us to compare the effects of curvature on the film boiling process. 


Consider the models of film boiling from the horizontal flat plate, horizontal cylinder, 
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and sphere as shown in Figure 6. As seen in Figure 6 T a thin vapor gap thickness 
6 and vapor dome are common to all the geometries. In all cuses, the liquid evap- 
orates and flows into domes which are assumed to form in a temporal periodic man- 
ner at regular Intervals \ or $*\ One coordinate is eliminated through geometric 
symmetry thereby enabling one to employ the stream function in solving ‘he momen- 
tum equations. 

Basic Assumptions 

For all geometries, the creeping flow assumptions aro made, the gap thickness is 
taken as constant, and surface capillary waves, while influencing the boundary, do 
not significantly alter the heat transfer or the laminar nature of the flow. In con- 
sidering the energy transport, radiation is neglected, the properties are evaluated 
at the film temperature = (T w + T s )/2), and convective terms are neglected in 
comparison with the conductive terms in the energy equation. The convective term in 
the last assumption are usually accounted for by replacing the latent heat of vapori- 
zation with a modified heat of vaporization of the form 

\* = \ x (l +£c p ATAi) (7) 

Justification for these assumptions is given in references 2, 5, and 7. 

Governing Equations and Boundary Condition 

With the previous assumptions the governing momentum and energy equations can be 
written as 


Efy * 0; V 1 2 * T = 0 (8) 

The operators E 4 and V 2 are given in references 2, 5, and 7 for each of the three 
coordinate systems under consideration, that is, for the flat plate, cylinder and 
sphere geometries. Reference 2 also presents a solution of equation (8) with a con- 
vective term left in the equation, resulting in a modified latent heat of vaporization, 
similar to equation (7). 

For the boundary conditions, the tangential vapor velocity to the interface is assumed 
to be zero, the interface to be at constant temperature T s , and the wall at a 
constant temperature T w . Also, uniform vaporization is assumed at the liquid- 
vapor interface and, of course, zero velocity at the heater surface. 

Constraints 

The above assumptions and boundary conditions do not lead to a tractable problem be- 
cause the vapor gap thickness, interlace velocity, and vapor dome spacing are un- 
knowns^— Thus, three additional constraints are required. 

(1) Static Force Balance Constraint. - The average pressure P fl on the vapor side 

of the liquid vapor interface results from the weight of the supported liquid, ambient 

pressure, and surface tension effects. This average pressure, P 8 , must be balanced 

by the vapor pressure field arising from the flow of vapdr into the vapor dome 



( 9 ) 


/ P, dA = P g A 
* A 

whore Pj is the pressure evaluated at the liquid-vapor interface. 

From this balance arises the two similarity parameters t and B 0 , discussed in 
the body of the report. 

(2) Interface Energy Balance Constraint . - The energy conducted across the vapor 
gap vaporizes the liquid and a balance must be maintained between the advancing 
evaporating interface and the incident energy. The energy balance may be written as 

-pXu 6 = q = -k 8T/8n |j (10) 

where n represents the coordinate normal to the interface. In reference 2, the bal- 
ance is augmented to include the effects of subcooling and radiation. 

(3) Optimization Constraint . - In determining the vapor dome spacing parameters, 
the energy transfer or rate of entropy production is assumed to be a maximum 

8h w /8L = 0} 8 2 h w /8L 2 < 0 (11) 

where L represents a characteristic wavelength or dome-spacing parameter and 
h w is the area weighted heat transfer coefficient: 

h\v = kfb^ft/ ^ t) + h D< A D/ A T> (12) 

The hypothesis asserts that the characteristic wavelength adjusts the size and spac- 
ing of the vapor dome to optimize the rate of heat removal from the surface. The 
wavelengths predicted by tills assumption compare favorably with the experimental 
data. 


With these three constraints, the system is now determinant. The analytical solu- 
tions are listed in table I and shown in Figure 3 along with tie -cueral solution; 
Table I also shows some special limiting cases. 

Heat Transfer Similarities 

The solutions for the horizontal flat plate (Ref. 2), cylinder (Ref, 5), and the sphere 
(Ref. 7) are presented in table I in terms of Rayleigh numbers defined as, 


R a d = P(P L “ P) d3 /(ty (13) 

Hj £ = p(p h - p)gX* 1 3 /(hp AT) (14) 


The flat plate represents the limiting case of a surface with no curvature, while the 
horizontal cylinder has one finite curvature. Except for the slight difference in the 
constants of proportionality (0. 424 compared with 0.373) the finite curvature of the 
wire increases the heat transfer by a factor 




1/4 


( 15 ) 






For small cylinders (small Bond numbers, B 0 ) the heat transfer coefficient will be 
significantly affected by the curvature of the surface, likewise for the sphere. Fig-* 
ure 3 illustrates this quite graphically. At small Bond numbers, figures 3 and 4 
show that the heat transfer coefficient from the sphere will be greater than that from 
an equivalent cylinder because the second-component of curvature on the sphere fur- 
ther enhances-ihe heat transfer. 


NOMENCLATURE 



c 

d 

a 

G(B 0 ) 

g 

Sc 


h 

k 


i 

L 

Nu 

P 

q 



AT 


area 

Bond number (d/t )2 
specific heat at constant pres- 
sure 

specific heat 
diameter 

derivative with respect to time 
differential operator 
function defined in reference 7 
local acceleration 
gravitational conversion con- 
stant 

heat transfer coefficient - — 
thermal conductivity 
effective thermoco uple the rmal 
conductivity, |fo/k^ +V k W 2 

dom e spacing para meter , 

»fiiSc/<PL- P)g 
characterizing length 
Nusselt number 
pressure 
heat flux 
radius 

radius of sphere 
modified Rayleigh number (Eq (13)) 
time 

quenching time 
temperature 

temperature difference (T w - T g ) 


u,v velocity components 

V volume 

6 vapor gap thickness 

p density 

a surface tension 


0 

0* 

X* 

h 

*c 

X 


stream function 


reference angle (Fig. 6) - values 
tabulated in ref. 7 


modified latent heal of vaporiza- 
, tioh (eq. (ll)) 
latent heat 


characteristic waveleng th 
wavelength 


t c system-time constant 
H viscosity 

Subscripts; 
d diameter 

D vapor dome 

f film conditions 


fb film boiling 


I interface 

l based on characteristic length 

L liquid 

Leid Leidenfrost 
s saturation 


T total 

w surface or wall 

w. Q initial surface or wall 
wire (thermocouple) 

6 evaporating interface 

y average, uninstrumented 










Table I. - Theoretical Heat Transfer Expressions and Limiting Cases for 
Three Geometries, The -Elat Plate, Cylinder, and Sphere 


Nusselt number expression 


Horizontal Flat Plate 
Nu { - 0.424 (Raj) 1 / 4 


Horizontal Cylinder 


kt hd A 

Nu h = -— = 0. 


Limiting cases 
A special limit in itself 


Large B 0 


0. 37 3,r„ 3 ^ 1/4 NUd . # ^l/« 


Sphere 


Small B 0 

/Rad\ 1/4 

Nu d = 0.00846^-^j 
. Large B_ 


Nu d = 2 + 


1 2Ra g — — j + 0 . mfRaJj^/S^A Nu d = 0. 35(Ra*j ^/ b“) 1/4 


4 CSC0*| (1 + COS-0 ) 


Small 


Nu d = 3 + (Ra3^) J 


.♦./tTV/4 


xpL2L + 0.177] 

by* J 



Theoretical comparison of curvature effects on pool film boiling for the 1/ 

cylinder, and sphere. The classic Bromley equation is for reference. 

Figured - Film boiling heal transfer data for spheres. 











